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ABSTRACT 

We present templates for the Sunyaev-Zel'dovich (SZ) angular power spectrum based on four mod- 
els for the nonlinear gas distribution. The frequency-dependent SZ temperature fluctuations, with 
thermal (TSZ) and kinetic (KSZ) contributions, are calculated by tracing through a dark matter sim- 
ulation, processed to include gas in dark matter halos and in the filamentary intergalactic medium. 
Different halo gas models are compared to study how star formation, energetic feedback, and nonther- 
mal pressure support influence the angular power spectrum. The standard model has been calibrated 
to reproduce the stellar and gas fractions and X-ray scaling relations measured from low redshift 
clusters and groups. The other models illustrate the current theoretical and empirical uncertainties 
relating to properties of the intracluster medium. Relative to the standard model, their angular power 
spectra differ by approximately ±50% (TSZ), ±20% (KSZ), and ±40% (SZ at 148 GHz) for I = 3000, 
ag = 0.8, and homogeneous reionization at z = 10. The angular power spectrum decreases in am- 
plitude as gas mass and binding energy is removed through star formation, and as gas is pushed out 
to larger radii by energetic feedback. With nonthermal pressure support, less pressure is required to 
maintain hydrostatic equilibrium, thus reducing the thermal contribution to the SZ power. We also 
calculate the SZ templates as a function of as and quantify this dependence. Assuming Ci oc (cr8/0.8)", 
the effective scahng index ranges from 7 < axsz ^ 9, 4.5 < axsz £ 5.5, and 6.5 < asz(148 GHz) < 8 
at / = 3000 for 0.6 < cts < 1- The template spectra are publicly available and can be used when 
fitting for the SZ contribution to the cosmic microwave background on arcminute scales. 
Subject headings: cosmology: theory - cosmic microwave background - large-scale-structure of universe 
- galaxies: clusters: general - intergalactic medium - methods: numerical 
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1. INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect imprinted in maps 
of the cosmic microwave background (CMB) is a promis- 
ing probe of the evolution of large-scale-structure. CMB 
photons propagating through the expanding universe 
are scattered by energetic electrons in the intraclus- 
ter medium (ICM) and intergalactic medium (IGM), 
resulting in secondary distortions which are the dom- 
inant temperature anisotropies on arcminute scales 
(|Sunvaev fc Zeldovichl[T970lll972| ). Measurements of the 
frequency-dependent distortions can be used to study 
the gas distribution in galaxy clusters and groups, and 
thus to understand the growth of structure. The electron 
scattering also traces the diffuse baryons in the filamen- 
tary cosmic web, which have so far been elusive in de- 
tection. Furthermore, the SZ effect has the potential to 
probe the epoch of reionization when the majority of elec- 
trons were dissociated from hydrogen a nd helium atoms. 
For reviews, see Birkinshaw (1999 ) and iCarlstrom et al.l 
§002). 

CMB experiments now have the sensitively and reso- 
lution to observe the temperature anisotropies with of 
order 10/iK noise and arcminute beams. The major 
SZ science includes the autocorrelation of temperature 
fluctuations, cross-correlation with large-scale-structure 
from galaxy surveys, and direct detention of galaxy clus- 
ters. In this paper we focus on the SZ angular power 
spectrum, an interesting statistic because its amplitude 
depends strongly on the normalization of matter per- 
turbations, generally parameterized by as- Assuming 



the SZ power scales as C'l oc Ug, previous theoreti- 
cal calculations find ax sz > 7 fe.g. ISeliak erani200l 
'Kom atsu fc Seliak|[200^ f or the t hermal SZ (TSZ) com- 
ponent, and axsz ^ 4 fe.g. lVishnia c 19871 for the kinetic 
SZ (KSZ) component. In principle, even an uncertain 
measurement only within a factor of 2 results in a better 
than 15% determination of as- 
Several groups have looked for the SZ contribution to 
the CMB angular power spectrum beyond the damp- 
ing tail at multipoles / > 1000. Experiments such 
as the Arcminute Cosmology Bolometer Array Re- 
ceiver (ACBAREI), Berkeley- Illinois-Maryland Associa- 
tion (BIMA0), and Cosmic Background Imager (CB0) 
report a significant detection of excess power coming 
from the SZ effect and extragalactic point sources. For 
example. iSievers et al.l p009il find an excess that is 1.6cr 
above the level expected for as = 0.8 from CBI observa- 
tions at 30 GHz. However, other experiments such as the 
Atacama Pathfinder Experiment (APEX-S^, Caltech 
Submillimeter Observatory (BolocanQ), Quest at DASI 
(QUaE0) and Sunyaev-Zel'dovich Array (SZ^H) report 
no large excess power. 
Most recently, the Atacama Cosmology Telescope 

^ http:/ /cosmology.berkeley.edu/group/swlh/acbar/ 

^ http://bima.astro.umd.edu/ 

^ http:/ /www. astro. caltech.edu/~tjp/CBI/ 

http:/ /bolo. berkeley.edu/apexsz/ 
^ http:/ /www. cso.caltech.edu/bolocam/ 
^ http:/ /en. wikipedia.org/wiki/QUaD 
^ http://astro.uchicago.edu/sza/ 
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(ACT0) and South Pole Telescope (SPTQ) have made 
unprecedented signal-to-noise measurements of the CMB 
angular power spectrum out to I ~ 10'*. Both groups 
fit for the SZ contribution using templates constructed 
for a ACDM cosm ology with erg — 0.8, which we pre- 
sented in Schg al et al.l ()2010f ) . They allow the normal- 
ization to vary and measure the s caling factor ^sz = 
Ci{crs)/Ci{0.8). iFowler et al.l (l2010l) place an upper limit 
(95% confidence level) of Asz < 1-63 based on ACT ob- 
servations at 148 GHz. Th is implies as < 0-8 6 (95% 
CL) if Ci oc (j\ is assumed. iLueker et al.l (|2009D report 
a best-fit value of A^z = 0.42 ± 0.21 (at 153 GHz) for 
SPT observations near 150 and 220 GHz. When com- 
bined with the Wilkinson Microwave Anisotropy Probe 
(WMAFH) 5-year constraints, the joint observations 
yield 0-8 = 0.773 ± 0.025. 

Currently, there are two publicly available sets of tem- 
plates for the f requency-dependent SZ a ngular power 
spectrum. The iKomatsu fc Seljakl (|200a KS02) tem- 
plateF^ is a popular choice and has been used in most 
analyses to date. The angular power spectrum is 
calculated using an analytical halo model, where the 
gas has a polytropic equatio n of state and is in hy- 
drosta tic equilibrium with a iNavarro. Frenk. fc Whira 
(jl997l NFW) gravitational potential. With thermal but 
not kinetic contributions to the temperature fluctua- 
tions, this template is expec ted to scale approximately as 
Ci oc al{Vt^,hf. In Sc heal el^ (|20Tnl) . we constructed 
SZ mapf[3 for an octant of the sky by tracing through 
a dark matter simulation processed to include gas in 
dark matter halos and in the filamentary IGM. The halo 
gas distribution is modeled with a polytropic equation 
of state, and the hydrostatic balance is performed di- 
rectly on the simulated dark matter halos to preserve the 
conc entration, substructure, and triaxiality of e ach sys- 
tem (jOstriker et al.|[2005l: iBode et al.|[2007l[2009l) . In the 
standard model o f IBode et al.l(|2009D . star formation and 
energetic feedback are calibrated against observations of 
low redshift clusters and groups. These templates have 
recently been used by ACT and SPT, as discussed above. 

In this paper, we address two main questions necessary 
to interpret observations. How do the SZ temperature 
fluctuations depend on the assumed model for the nonlin- 
ear gas distribution? How does C; scale with as for each 
given model? We present templates for the SZ angular 
power spectrum based on four models for the nonlinear 
gas distribution in a ACDM universe with arbitrary cg. 
Different halo gas models are compared to study how star 
formation, energetic feedback, and nonthermal pressure 
support influence the SZ effect. Section [2] reviews the for- 
malism for calculating the SZ temperature fluctuations, 
and Section [3] describes the construction of the numeri- 
cal models. The SZ angular power spectra are compared 
in Section 131 and constraints on erg from recent observa- 
tions are discussed in Section jS] We adopt the cosmolog- 
ical parameters: n,-a = 0.264, I^a = 0.736, rJb = 0.044, 
h — 0.71, Us = 0.96, and erg = 0.8 0, which are siniilar t o 
the recent WMAP 7-year results (|Komatsu et al.ll2010D . 

* http:/ /www. physics. princeton.edu/act/ 

^ http://spt.uchicago.edu/spt/ 

http:/ /lambda. gsfc.nasa.gov/product/wmap/ 

http:/ /lambda. gsfc.nasa.gov/product/map/dr4/pow_sz_spec_g( 

http:/ /lambda. gsfc.nasa.gov/toolbox/tb_cmbsim_ov.cfm 



2. SUNYAEV-ZEL'DOVICH EFFECT 

The SZ effe ct is commonly considered to have tw o main 
components ( Sunvaev fc ZeldovichI Il97a I1972D . The 
thermal SZ (TSZ) term arises from inverse Compton 
scattering of the CMB with hot electrons, predominantly 
associated with shockheated gas in galaxy clusters and 
groups. The kinetic SZ (KSZ) term is a Doppler term 
coming from scattering with electrons having fast, pe- 
culiar motions. Another common distinction is that the 
KSZ effect has a nonlinear component associated with 
the small-scale ICM, and a more linear component com- 
ing from the large-scale IGM. The signal arising from 
the latter was flrst calculated for the linear regime by 
lOstiikeiifc Vishniac (1986) and Vislmiac (1987), and is 
often referred to as the OV effect. In this section, wc flrst 
write down the formalism for the nonrelativistic limit and 
then the more general relativistic case. 

Modeling the SZ effect requires knowing the number 
density rie, temperature Tg, and velocity Ve of the elec- 
tron distribution. In the nonrelativistic limit, the change 
in the CMB temperature at frequency v in the direction 
h on the sky is given by 

^„,, = (^) J^) ,1, 

where the dimensionless Compton y and Doppler b pa- 
rameters, 



UeTedl = / 6'edT, 



6 = ^ / UeViosdl = / Aosdr, 



(2) 



(3) 



are proportional to integrals of the electron pressure 
and momentum along the line-of-sight (los) , respectively. 
The dimensionless temperature, los peculiar velocity, and 
the optical depth through a path length dl are given by 



1.96 X 10" 



— = 3.34 X 10-^ 
c 



keV 



100 km/s 



(4) 



(5) 



dr = 



arn^dl = 2.05 x 10"^ f ) ( JL] , (6) 

VlO-3 cm3/ ^Mpc J 

respectively. For a typical cluster, the Compton y param- 
eter is expected to be approximately an order of magni- 
tude larger than the Doppler b parameter. 

The nonrelativistic TSZ component has a frequency 
dependence as specifled by the function. 



= x„ coth(a;i./2) 



= hv/{kBTcMB), (7) 



which has a null at « 218 GHz. The distortion ap- 
pears as a temperature decrement at lower frequencies 
and as an increment at higher frequencies relative to the 
null. The nonrelativistic KSZ component is independent 
of frequency, but the sign of the distortion depends on 
the sign of the los velocity. We chose the convention 
..cfnwhere vios > if the electrons are moving away from the 
observer. 
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In the general, relativistic case, the change in the CMB 
temperature at frequency v in the direction n is given by 



AT 



T, 



CMB 



(n) 



-Yo + -Fi 
6 " 3 



dr. 



(8) 



where the y's and C"s are know n frequency-dependent 
coefficients (N ozawa et alJ lTQQSf). Note that Eq. [T] con- 
tains only the first-order terms in Eq.[S]and that fi, — Yq. 

In this paper, we use Equation ([8]) to calculate the total 
SZ signal. When the TSZ and KSZ effects are discussed 
individually. Equations ([2]) and are used to calculate 
the Compton y and Doppler 6, respectively. After multi- 
plying the Compton y by Tcmb = 2.726 /xK, it becomes 
equivalent to the nonrelativistic [ATtszI at 146 and 280 
GHz. The Doppler b temperature fiuctuation is equiva- 
lent to ATksz at all frequencies. 

3. SIMULATIONS 

The SZ effect is modeled by postprocessing a simula- 
tion of the large-scale structure of the universe, as de- 
scribed in detail in Se hgal et al.l (|201Clf ) . In this section, 
we summarize the methodology and present additional 
details on the modeling. 

The matter distribution along a past light cone, span- 
ning the redshift range < z < 10 and covering one 
octant of the sky, is generated from an N-body simula- 
tion of 1024^^ dark matter particles evolved within a pe- 
riodic box of comoving side length L = 1000 h~^Mpc. 
The resolution is set by the particle mass, nip = 
6.82 X 10^" h~^MQ, and the gravitational spline soft- 
ening length, e — 16.28 h'^kpc. The friends-of-friends 
(FoF) halo mass functio n is in good agreement with 
the fitting formula from iJenkins et al.l ([200 Ij down to 
AfpoF 7 X 10^2 (100 particles). 

For lower redshifts z < 3, the positions and velocities of 
all particles in the light cone are saved for postprocessing. 
The particles thus saved allow a higher-resolution recon- 
struction of dark matter halos and filaments in the large- 
scale structure. At higher redshifts 3 < z < 10, projected 
information is saved instead. Within thin redshift shells, 
particles are subdivided by angular coordinates and then 
projected along the line of sight to construct surface den- 
sity fields for mass and momentum. There are 579 such 
shells, with thickness Az w 0.09 at z w 10 and Az w 0.03 
at z ~ 3. 

The corresponding gas and electron distributions are 
modeled with contributions from three components: 

1. Gas in massive dark matter halos with MpoF > 
2 X 10-^^ h~^MQ at z < 3 is modeled with a poly- 
tropic equation of state an d in hydr ostatic equilib- 
rium (.Ostriker et al. 2005 : .Bode et a l. 2007, 2009). 
For each halo, the matter profile is reconstructed 
at high resolution using the saved N-body parti- 
cles. This procedure preserves the concentration, 
substructure, and triaxiality of the system. See 
Section [3?T] for additional discussion. 



2. Gas in lower mass halos and the IGM at z < 3 are 
modeled using all other saved N-body particles not 
associated with the massive halos discussed above. 
The effective gas velocity and temperature are ap- 
proximated using the peculiar velocity and velocity 
dispersion of the particles, respectively. 

3. At higher redshifts 3 < z < 10, the gas is as- 
sumed to trace the matter, and the temperature 
(T ~ 10'' K) is predominantly set by photoioniza- 
tion rather than shockheating. The electron distri- 
bution is modeled using the saved projected mass 
and momentum density fields. Hydrogen rcioniza- 
tion is assumed to occur instantaneously at z = 10 
and helium is only singly ionized. 

Sky maps of the SZ effect are made by tracing through 
the simulated electron distribution and projecting the ac- 
cumulated tem perature fluctuations onto a HEALPi30 
(|G6rski et al.l |2005) grid with pixel resolut i on of 0.4 
arminutes (Nside = 8192). See ISehgal et al.l (|2010l) for 
additional details and discussion. 

3.1. Halo gas models 

The hot gas distribution associated with galaxy clus- 
ters and groups is modeled as having a polytropic 
equation of state and being in hydrostatic equilibrium 
(jOstriker et al.|[2005l: [Bode et al.l !2007. 2009). In this pa- 
per, we compare four gas models to study how gas physics 
influence the SZ effect, in particular the angular power 
spectrum of the secondary temperature anisotropies: 

1. Adiabatic model, which has neither star formation 
nor energetic feedback to represent nonradiative 
gas physics. 

2. Standard model, with star formation and feedback 
calibrated against observations of nearby clusters 
and groups. 

3. Low- /gas model, with twice the star formation and 
twice the feedback energy input as the standard 
model. 

4. Nonthermal20 model, which includes 20% nonther- 
mal pressure support plus more star formation, 
but less feedback energy input, than the standard 
model. 

Figure [T] shows the stellar mass fraction and gas 
mass fraction /gas within i?5oo at z < 0.2 for these mod- 
els. The mass fractions are defined relative to M^qq and 
normalized by the cosmic baryon fraction fib/^^m- The 
adiabatic model provides an appropriate upper limit on 
the gas fraction. In the absence of star formation and 
feedback, the gas fractions are close to the cosmic aver- 
age and varies only weakly with mass and redshift. The 
slight increase towards lower masses is due to the higher 
halo concentrations. While this basic model is known to 
be incorrect, it is still useful in comparison for undcr- 
st anding the effects of radiative gas physics. 

IBode et"aLl (2009.) previously constructed the stan- 
dard model to match the stellar and gas fractions 

^'^ http://healpix.jpl.nasa.gov/ 
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Fig. 1. — Top: Gas fractions for halos at 2 < 0.2 from the stan- 
dard (black solid), adiabatic (red long-dash), low- /gas (blue short- 
dash), and nonthermal20 (green dotted) models. The Icr (dark 
gray) and 2(t (light gray) scatters from the stan dard model are 
also shown. Obser vational data po ints are from IVikhlinin et al.l 
([2006, circles) and I Sun et al.1 ||2009| . triangles). Bottom: Stellar 
fractions for the same halo s as above . The standard model adopts 
the best-fit relation from I Lin et aP {2003, circles), the low- /gas 
model has twice the stellar fraction as the standard mode l, while 
the nonthermal20 model is taken from lGiodini et al.l 1)20091 ). 



observed in nearby clusters and groups. The stel- 
lar fractions at z — are g iven by f^,{< R500) = 
0.0164(Af5oo/3xl0i4 Mq)-°-'^'^ (|Lin et al.ll2003D . and the 
redshift evolution follows a delayed exponential model 
(jNagamine et al.l [2006). In ord er to match the gas 
fractions from X-ray observations (jVikhlinin et al.lf2006l : 
iSun et al.| [2009'). a feedback efficiency e = 4 x 10"^ and 
energy input of eMpc^ is required, where Mp is the mass 
of formed stars. With this calibration, the standard 
model reproduces several observed X-ray scaling rela- 
tions, such as Tx - M500, Lx - M500, and Yx - M^qq. 
We choose the standard model as the fiducial case when 
comparing with other models and with observations. 

The low- /gas model has even higher stellar fractions 
and lower gas fractions than the standard model. This 
model has twice as much star formation and twice as 
much feedback energy input. The feedback efficiency 
constant is the same for both models, and for the adopted 
value we find that star formation and feedback have 
comparable effects on reducing the gas fraction within 
i?5oo- The gas fractions skirt the lower bounds of the 
uncertainty in the measurements from nearby clusters 
and groups. While the Tx — M500 relation for M500 ^ 
lO^'' Mq is in good agreement with observations, the Lx 
and Yx values are lower by approximately 20%. The 
low- /gas model will help in understanding what is a rea- 
sonable lower bound on the SZ angular power spectrum 
due to uncertainty in the star formation and feedback 
histories. 



The nonthermal20 model is different from the other 
three models in that 20% of the hydrostatic pressure 
is assumed to be nonthermal. Turbulent and relativis- 
tic contributions account for 15% and 5%, respectively, 
of the pressure at all radii. Note these two compo- 
nents contribute differently to the total energy budget, 
because of the different ratios of specific heats. This 
model also includes star formation and feedback, but 
different in detail than previously discussed. The ste l- 
lar fraction follows Equation 10 of iGiodini et al.l (|2009t ): 
/,(< R500) = 0.05(M5oo/5 X 1013 Mo)-o-26. This is a 
fit to X-ray-se lected groups an d poor clusters of lower 
mass than the iLin et al.l (|2003[ ) sample; as it turns out, 
this is almost twice as large as the /, used in the stan- 
dard model. A lower feedback efficiency, e = 1 x 10~^, 
is needed to reproduce the gas fractions and scaling rela- 
tions from X-ray observations. The model does this well, 
keeping in mind that the total masses derived observa- 
tionally must also be corrected for the 20% nonthermal 
contribution in the hydrostatic equilibrium equation. 

There is observational support for including both tur- 
bulent/bulk motions and cosn iic rays. Co nrparing X- 
ray and weak lensing masses, iZhang et al.l ( |2010) find 
nonthermal pressure support of ~ 9% inside -R500, while 
iMahdavi et al.l (|2008f ) found roughly twice this amount. 
Comparing the gravitational potential profiles of the 
central galaxies of the Fornax and Vir go clusters (de- 
rived from X-ray and from optical datal. lChurazov et al.l 
(jlOOS) find the total nonthermal pressure to be < 10 — 
20% of the gas thermal pressure. Examining the width 
of X-ray emission lines. I Sanders et aD (|2010f ) place an 
upper limit of 13% on the turbulent energy density as 
a fraction of the thermal energy density in the core of 
A1835. Simulations of clusters give values for energy 
in cosmic rays relati ve to the thermal en ergy of approxi- 
mately 5 — 10% fe.g. lJubelgas et aLll2008[ ). in good agree- 
ment with con straints from gam ma ray and radio ob- 
servations (e.g. iPfrommer fc Enfi lin 2004). Simulations 
likewise give contributions from bulk flows and turbu- 
lence of approximately 5 — 20% (e.g. iLau et al.l 120091 : 
iMeneghetti et al.ll2009l : iBurns et al.ll2010D . ThuTto ex- 
amine the effects of these components we include a gen- 
erous 20% nonthermal pressure support in the nonther- 
mal20 model to examine the effects of these components, 
with turbulence being the dominant contributor. 

Figure [2] shows the radial pressure profiles for halos 
with M500 > 10^^ Mq and z < 0.2 drawn from our mod- 
els, and compares them with other theoretical and ob- 
servational results. We rescale t he halo profiles b y the 
characteristic pressure P500 (see iNagai et al.l 120071 ) and 
calculate an average by weighting each h alo by th e inte- 
grated pressure within i?5oo (see Battag lia et al.| [2010'). 
At smaller radii {R < 1.5i?5oo), the differences in our 
models are due to the depletion of gas mass and binding 
energy through star formation, the redistribution of gas 
to larger radii through energetic feedback, and the reduc- 
tion in the amount of thermal pressure needed to main- 
tain hydrostatic equilibrium in the presence of a non- 
thermal component. At larger radii (i? > 1.5i?5oo), ha- 
los experiencing more energetic feedback have shallower 
profiles as gas gets pushed outwards. When consider- 
ing the total pressure for R < i?vir ~ 2i?5oo, the dif- 
ferences between the standard, adiabatic, and low- /gas 
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Fig. 2. — Left: Average pressure profiles for hales with M^qq > 10^* Mq and z < 0.2 from the standard model (black solid), adiabatic (red 
long-dash), low- /gas (blue short-dash), and nonthermal20 (green dotted) models. The la (dark gray) and 2(t (light gray) scatters from the 
standard model are also shown. The bottom panel shows pressure ratios with respect to the stand ard model. No t e a ki nk at R/R500 ~ 3 
is present because the IGM contribution is not included here for the simulated profiles. Right: The lArnaud et all 120091 . cyan dotted) and 
IKomatsu &: SeliakI l(200ll . [2003. orange short-long-dash) profiles are for M500 = lO^"* Mq (lower) and lO^'' (upper) at z = 0.2. The 
IBattaglia et al.l 1)20101 . magenta dot-dash) profile is a weighted average from simulated halos in the mass range 10" < Msoo/Mq < 1015 at 

2 = 0. 

models are almost entirely due to the different stellar 
fractions. In the case of the nonthermal20 model, there 
is an additional reduction due to the 20% nonthermal 
hydrostatic support. The pressure profiles rapidly de- 
crease at i?/i?5oo ^ 3, as the IGM contribution is not 
included in the plot. The gas in each halo has a differ- 
ent, anisot ropic outer limit, w hich also changes with the 
model fseelSehgal et al1 l20lQV 

lArnaud et al.l (120091) used the REXCESS cluster sam- 
ple (jBohringer et al'fl2007[ ) at z < 0.2 and found that the 
scaled pressure profiles are well fit by a generalized NEW 
pressure model pagai et al.|[2007l ). The pressure profiles 
from the standard model are in good agreement with this 
recent analysis. At small radii {R/Rsoa < 0.1), we under- 
predict the pressure because of the finite resolution of the 
gas reconstruction, which is computed using a Cartesian 
grid with cell length Ax = 32.55 /i~^kpc (comoving) « 
0.03i?5oo- However, this small region contributes only 
minorly to the total pressure. At intermediate radii 
(0.1 < R/R500 < 1), which is the most relevant range 
for this comparison, our profiles ar e similar to the REX - 
CESS sample (also see Eigure 8 in lArnaud et all [20091 ) . 
Note that on larger radii (1 < R/R500 < 4), their best-fit 
model is based on r esults from hydrodynamic simulations 
(iNaeai et al.ll2007D rather t han X-ray d ata. 

The IKomatsti fc SeliakI (pfM \2(m . KS02) pressure 
profiles shown in Eigure [2] are for M500 = 10^** Mq and 
10^5 Mq at z = 0.2 (courtesy of E. Komatsu). While 
their gas model is also based on a polytropic gas in hy- 
drostatic equilibrium, there are some important differ- 
ences to note. Their model is based on nonradiative gas 
physics (having no stars nor feedback), but the pressure 
profiles differ from our adiabatic model for the follow- 



ing reasons. Eirst, the KS02 model is normalized such 
that the baryon fraction at the virial radius is equal to 
the cosmic average, which has the effect of setting the 
baryon fraction within this radius to be less than the cos- 
mic average. This is more similar to the standard model, 
rather than the adiabatic. Second, the polytropic index 
depends on the halo concentration and is approximately 
1 . 1 for the mass (M500 > IO^-^Mq) and redshift {z < 0.2) 
ranges being considered in Eigure [U which is lower than 
our constant value of 1.2. As a result, the KS02 pressure 
profiles for clusters are shallower and do not decrease as 
fast with r adii. 

Battagl ia et al.l (|2010( ) recently modeled the SZ effect 
using hydrodynamic simulations with radiative cooling, 
star formation, and feedback from supernovae and ac- 
tive galactic nuclei ( AGN) . The simulations directly cap- 
ture bulk and turbulent gas motions, an d also i nclud e 
a prescription for cosmic rays (Pfromm er et al.l l2007| ) . 
In Eigure [21 their scaled pressure profile (courtesy of 
N. Battaglia) is a weighted average from halos in the 
mass range lO^^ < Mgoo/Mg < lO^^ at z = 0. We find 
very good agreement with our standard model, within 
~ 20% (which is smaller than the la scatter) over the ra- 
dial range 0.1 < R/R500 < 3. The nonthermal20 model, 
which also accounts for cosmic rays and turbulent gas 
motions, compares even more favorably with theirs at 
R/R500 > Ij but less well at smaller radii. 

4. RESULTS 

We make sky maps of the SZ effect in Healpix for- 
mat with pixel resolution of 0.4 arcminutes (Nside = 
8192). Eigure [3] shows a 4 degree x 4 degree sample field, 
centered on a cluster with A/vir ~ lO^^M^ a,t z ^ 0.3. 
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Fig. 3. — Top left: The SZ temperature anisotropies in fiK at 148 GHz from a 4° X 4° sample field in the standard model map. At this 
frequency, the TSZ signal is a temperature decrement, but the KSZ signal is a decrement (increment) if the los peculiar velocity is positive 
(negative). The cluster at the center of the field has mass M^ii ~ 10^'^ Mq and redshift z Ri 0.3. Top right: The temperature difference, 
^Tlidiabatic ~ ^T'st^ndard i between the adiabatic and standard models. In the adiabatic model, the gas is more centrally concentrated, 
leading to relatively stronger central decrements, but slightly weaker decrements at larger radii. Bottom left: The temperature difference, 
^Ti^^_f^^^ — ^Tljtandardi between the low- /gas and standard models. With more star formation and feedback, the low- /gas model has weaker 
central decrements, but slightly stronger decrements at larger radii. Bottom right: The temperature difference, Arnontjn,rmal20~^^standardi 
between the nonthermal20 and standard models. The nonthermal20 model has less pressure and weaker decrements at all radii. 
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Fig. 4. — Cumulative probability distribution functions of Comp- 
ton y (right) and Doppler b (left) temperature fluctuations for the 
standard (black solid), adiabatic (red long-dash), low- /gas (blue 
short-dash), and nonthermal20 (green dotted) models. The y and 
b values have been multiplied by Tcmb = 2.726 /iK to get the tem- 
perature fluctuations in fiK. The angular smoothing scale is set by 
the map pixel resolution of 0.4 arcminutes. 

The top-left panel shows the SZ temperature anisotropics 
in /j,K at 148 GHz from the standard model, while the 
other panels show the temperature differences between 
the standard model and the other three models for the 
same sample field. At this frequency, the TSZ compo- 
nent is always a temperature decrement, but the KSZ 
component can be positive or negative. 

For comparison, the reference cluster in the middle of 
the field has a central decrement of approximately 355, 
280, 225, and 200 /iK in the adiabatic, standard, low- /gas, 
and nonthermal20 models, respectively. Comparing the 
first three models, we find that the central decrement 
gets weaker as more gas is removed through star forma- 
tion and as more gas is pushed out to larger radii by 
energetic feedback. The latter also results in relatively 
stronger decrements at larger radii. Compared to the 
standard model, the larger stellar fraction and the 20% 
nonthermal pressure support in the nonthermal20 result 
in less pressure and weaker decrements at all radii. In 
Sections 14.11 - 14.31 we show that the effects of star for- 
mation, energetic feedback, and nonthermal hydrostatic 
support lead to scale-dependent changes in the power 
spectrum. 

Figure H] compares the cumulative probability distri- 
bution functions (PDF) of Compton y and Doppler b 
temperature anisotropics from maps with pixel resolu- 
tion of 0.4 arcminutes. After multiplying the values by 
Tcmb = 2.726 /xK, the Compton y becomes equivalent to 
IATtszI at 146 and 280 GHz, and the Doppler b equiv- 
alent to ATksz at all frequencies. On average, the \b\ 
temperature fluctuations are approximately an order of 
magnitude smaller than y. In principle, the PDF of the 
SZ signal can be used in place of the halo mass function 
to constrain cosmological models. However, it will be 




Fig. 5.— Top: Compton y or TSZ (at 146 and 280 GHz) angular 
power spectra for the standard (black solid), adiabatic (red long- 
dash), low- /gas (blue short-dash), and nonthcrmal20 (green dot- 
ted) models. For comparison, the Komatsu & Seljak (2002, orange 
short-long-dashes) and Battaglia ct al. (20f0, magenta dash-dot) 
templates have been converted to the same frequency. Bottom: 
Ratios of power spectra with respect to the standard model. Al- 
though the KS02 template is based on nonradiative gas physics, 
it resembles the standard model much more than the adiabatic 
model. 

difficult to break the degeneracy between cosmological 
and astrophysical parameters, for example between 
and /gas(Af, z), using this statistic alone. Furthermore, 
it will be difficult to isolate the SZ one-point signal in 
the presence of the CMB, extragalactic point sources, 
galactic dust, and noise. 

4.1. TSZ angular power spectra 

Angular power spectra of the Compton y (and equiv- 
alently the TSZ at 146 and 280 GHz) temperature fluc- 
tuations are shown in Figure [H The power spectra are 
calculated from the simulated maps using the software 
PolSpiccH (e.g. ISzapudi et aH I200TI : IChon et al.l [200I 
and have been corrected for the pixel window. In all 
four models, the TSZ power spectrum is broadly peaked 
at I ^ 3000, where it also becomes dominant over the 
lensed CMB power spectrum. 

Relative to the standard model, the adiabatic model is 
47% larger, the low- /gas model is 30% smaller, and the 
nonthermal20 model is 45% smaller at / — 3000. The ef- 
fects of star formation and feedback can be simply under- 
stood. As more gas mass and binding energy is removed 
through star formation and as more gas is pushed out 
to larger radii by energetic feedback (see Figure [3]), the 
overall amplitude of the TSZ angular power spectrum de- 
creases. Since feedback more strongly affects lower mass 
halos because of their shallower potential wells, there is 
a stronger reduction in power on smaller angular scales. 
Correspondingly, as more massive or lower redshift halos 
make a relatively larger contribution to the total power, 

http:/ /www. planck.fr/articlef41. html 
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Fig. 6. — Left: The fraction of the thermal power (from both halos and the IGM) at I = 3000 as a function of minimum mass and 
maximum redshift in the standard model (black solid). The nine curves (from bottom to top) arc for PDF(> M500, < z) = [0.1, 0.9] in 
increments of 0.1. Right: PDF(> M500, < z) = 0.5 (bottom) and 0.9 (top) from the adiabatic (red long-dash), low- /gas (blue short-dash), 
nonthermal20 (green dotted), and KS02 (orange short- long-dashes) models are compared to the standard model. 



the peak in the spectrum shifts to larger angular scales. 

Interestingly, we find that the nonthernial20 model has 
even lower power than the low- /gas model at all angular 
scales, despite having higher gas fractions. Since their 
stellar fractions are very similar, the removal of gas mass 
and binding energy through star formation is very close 
to equal in these two models, and thus not responsible 
for the differences. The two reasons for the differences in 
the power spectra are as follows. First, while energetic 
feedback in the low- /gas model redistributes the mass 
and energy to larger radii, the nonthermal hydrostatic 
support in the nonthermal20 model simply reduces the 
pressure at all radii (see Figure [S]). Second, the feedback 
energy per baryon is correlated with the stellar fraction, 
which decreases at higher redshifts, while the assumed 
20% nonthermal support is independent of z. For reason- 
able values of parameters, nonthermal support can have 
comparable effects as energetic feedback on lo wering the 
ampli tude of the thermal power. Also, see iShaw et al.l 
(|2010( 1 for recent work using halo model calculations to 
study the impact of cluster physics on the TSZ angular 
power spectrum. 

Figure ini shows the fraction of the thermal power (from 
both halos and the IGM) at I = 3000 as a function 
of minimum mass and maximum redshift, denoted by 
PDF(> M500, < z). For the standard model, all halos 
with M500 > 1.5 X 1014 (3 X 10") Mq and z < 0.85 
(1.5) contribute approximately 50% (90%) of the power. 
The other three models require similar mass and redshift 
ranges, but there are systematic differences. Models hav- 
ing higher stellar fractions at low redshifts (e.g. low- /gas 
and nonthermal20 models) have relatively larger con- 
tributions coming from higher redshifts as /* declines. 
Also, models with gas fractions which decrease faster 
toward low masses (e.g. low- /gas) have relatively larger 



fractions of the power coming from massive halos with 
M500 ^ 3 X 10^4 Mq. Figure m demonstrates that in order 
to build templates for the TSZ angular power spectrum 
accurate to 10%, we need to calibrate theoretical pres- 
sure profiles using a fair sample of clusters and groups 
out to z ~ 1.5. 

Overall, the four models are s imilar to result s from 
hydrodynan iic sim ulations (e.g. iRefreeier et al.l 120001: 
Seliak et al.l 120011: iSprineel et al.l 120011: Ida Silya et alJ 



200H IBondeTali 120021 120051 iHallman et al.l I2007t 

Battaglia et al.| l2010')" and semi-analytical/numeric al caL 
culat i ons (e.g. Komatsu & Scljak 2 002; Holder et al 
2007 [ IMoodlev^t al.ll200ft iPeel et al.l[2rM IShaw et al 
2010D . However, there are still factor of 2 differences 



amongst the various predictions, even after correcting 
for different cosmologi es using the a pproximate scaling 
Ci (X al{n\,hf (Selia k et aiTl2001t Komatsu & ScljaSj 
|2002|) . It is generally difficult to correct for different as- 
trophysical assumptions relating to star formation and 
feedback. Note that previous hydrodynamic simulations 
are often limited to small box sizes, where the abundance 
of clusters (and even groups) are underpredicted due to 
missing large-scale power in the initial conditions. Thus, 
these results are highly prone to sample variance. 

For comparison, the popular KS02 template is also 
shown in Figure [5] It has been converted to the same fre- 
quency, but no correction for cosmological parameters is 
applied, because us = 0.8 is the same and {Vt^hY is larger 
by only 4% compared to our cosmology. At I = 3000, it 
has approximately 10% more power than the standard 
model. Although the KS02 model has no stars and feed- 
back, it resembles the standard model much more than 
the adiabatic model. This is due to the choice of den- 
sity normalization and polytropic indices, as previously 
explained in Section 13.11 The KS02 template also has 
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a slightly different shape, with relatively more power on 
larger angular scales for the following reasons. First, the 
massive halos have lower polytropic indices and more ex- 
tended pressure profiles. Second, when calculating the 
angular power spect rum with the Lim ber equation (see 
Equations (l)-(2) in lKomatsu &: Seliak 2002). the pres- 
sure profiles are integrated out to > 3i?vir to improve nu- 
merical convergence. Figure [6] shows that all halos with 
M500 > 5 X 10^3 X 1013) Mq and z < 1.1 (2.0) con- 
tribute approximately 50% (90%) of the thermal power 
at I = 3000 in the KS02 template. In their gas model, the 
concentration parameter and polytropic index decrease 
toward higher redshift at fixed mass, which enhances 
the pressure at R/R500 0.5. Consequently, high red- 
shift and low mass halos contribute relatively more to 
t he power than in any o f our models. 

iBattaglia et all ()2010f ) have calculated the TSZ angu- 
lar power spectra from many realizations of maps (max- 
imum size = 3.2° X 3.2°) made with hydrodynamic sim- 
ulations (maximum box size = 330 /i~^Mpc). In Fig- 
ure [5] we show their fiducial case based on a gas model 
with AGN feedback (courtesy of N. Battaglia). No 
correction for cosmological parameters is applied, since 
cTg = 0.8 is the same and (J7b^)^ is only 2% smaller 
than our value. At I = 3000, this case has approxi- 
mately 32% less power than the standard model, but 
23% more power than the nonthermal20 model. These 
differences are consistent with those found in the aver- 
age pressure profile of low redshift clusters (see Section 
13.11) . which is reassuring. However, it is not conclusive 
that the two approaches would yield similar power spec- 
tra for the same low redshift pressure model since the 
majority of the power comes from higher redshifts. Dif- 
ferent implementations of gas prescriptions may not have 
the same evolution of gas and stellar properties. Fur- 
thermore, the finite box and map sizes may have differ- 
ent amounts of sample variance, in particular at larger 
angular scales. A more detailed comparison, for exam- 
ple using the PDF(> M500, < z), is required in order to 
check the robustness of agreement in various approaches. 

4.2. KSZ angular power spectra 

Angular power spectra of the KSZ temperature fluc- 
tuations are shown in Figure [T] The KSZ signal is 
independent of frequency and its angular power spec- 
trum is much broader, with a peak amplitude that is 
approximately several times smaller compared to the 
Compton y angular power spectrum. This is consis- 
tent with results fro in hydro dy namic simulations (e.g. 



Springel et all 120011: |da Silva c t~all |200U iZhang et al.l 
20041: iHallma n et al.l \2009i. Battag ha et al.l \201&i a^ 



Hallma n et al.l [200^ Battag lia et al.l 
semi-analytical/numerica l ca lculati ons (e.g. [M a & Fry 



I20 

.Ju 



2002t [Zhang et al.l 120041: iHernandez-Monteagudo fc Hoi 
20091 ). Note that previously when erg ~ 1 was assumed, 
the KSZ term was expected to be over an order-of- 
magnitude smaller than the Compton y, but the ratio 
is closer with the currently favored value of cts 0.8. At 
I - 3000, we find Cksz/Ctsz ~ 0.19 (adiabatic), 0.23 
(standard), 0.29 (low-/gas), and 0.39 (nonthermal20), 
where the thermal term is evaluated at 146 or 280 GHz 
(i.e. Compton y). 

Compared to the standard model, the adiabatic model 
is 20% larger, the nonthermal20 model is 6% smaller, and 



1= 



1.5 



1 1 1 1 1 1 1 1 1 1 
^-^ 


1 1 1 

- 


/ 




/ 


^ 


It , ' 


- 


J ^ 




h' 




I 




— \ — \ — — \ — \ — \ — — \ — \ — \ — — ^ 


H 1 1 1 — 




y 

" , , 1 , , , 1 , , , 1 





2000 



4000 



6000 



8000 



Fig. 7. — Top: KSZ angular power spectra for the standard (black 
solid), adiabatic (red long-dash), low- /gas (blue short-dash), and 
nonthermal20 (green dotted) models, in which hydrogen is homo- 
geneously reionizcd at ^rcion = 10. For comparison, the standard 
model is modified to have hydrogen reionizcd later at Zrcion = 6 
(orange short-long-dash). Bottom: Ratios of power spectra with 
respect to the standard model. 

IGM contribution is the same for all of the models by 
construction, the differences in the power spectra come 
from the halo contribution alone, and since the models 
have the same velocity field the differences are entirely 
due to the halo gas fractions (see Figure [ij. Thus, the 
power spectrum from the adiabatic model has the highest 
amplitude while the low- /gas model has the lowest. The 
nonthermal20 model does not have the lowest amplitude 
for the KSZ effect, unlike for the TSZ effect. 

In our models, the KSZ signal includes scattering with 
the nonlinear electron distribution up to the epoch of 
cosmic reionization. In the fiducial case, hydrogen is as- 
sumed to be homogeneously reionizcd at ^reion = 10. The 
corresponding Thomson optical depth for electron scat- 
tering, tt = 0. 08, is consistent with the recent WMAP 
7-year results (iKomat su et al.l [20101 ). For the standard 
model, we also consider another case where hydrogen is 
reionizcd later, at Zreion = 6, to show the uncertainty 
due to the poo r ly constrained reionizatio n history (see 
iFan et al.l [20061 IMorales fc Wvith^ [20091 for recent re- 
views). Interestingly, we find that the uncertainty due 
to the halo gas fraction is comparable to the uncertainty 
due to the reionization history. Thus, in order to con- 
strain the reionization epoch through the KSZ angular 
power spectrum, it is clearly necessary to understand 
the nonlinear contribution from collapsed objects. The 
nonlinear KSZ and TSZ effects can be studied by cross- 
correlating maps of the CM B with those from galaxy 
redshift surveys (e.g. IHo et a l. 2009; Shao et al. 20l3). 

We have neglected to include the contribution from the 
inhomogeneous reionization epoch. This component is 
expected to be comparable or even larger th an the contri- 
bution from the fully reionizcd universe (e.g. Santos et al 



the low- /gas model is 10% smaller at I = 3000. Since the [2003; .McQuinn et al...2005; .Zahn et al...2005; .Iliev et al 
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Fig. 8. — Top: SZ angular power spectra at 148 GHz for the stan- 
dard (black solid), adiabatic (red long-dash), low- /gas (blue short- 
dash), and nonthermal20 (green dotted). The KS02 (orange short- 
long-dash) template has been scaled to 148 GHz and it does not in- 
clude any contribution from the KSZ effect. The AC T upper limit 
is derived using j4sz,STD < 1.63 (95% CL) taken from lFowler et aLl 
(I2OIOI) . The SPT best-fit value is based on Asz.STD = 0.42 ± 0.21 
taken from ILueker et al.l l|2009l) . This measurement was made at 
153 GHz, but the plotted power has been rescaled to 148 GHz. Bot- 
tom: Ratios of power spectra with respect to the standard model. 



120071) . It is especially important to account for this high 
redshift component when examining the angular power 
spectrum near the TSZ null frequency. However, at fre- 
quencies < 150 GHz and > 280 GHz, the KSZ contri- 
bution to the SZ angular power spectrum becomes less 
important as it is expected to be smaller than the current 
uncertainty in the TSZ term. 

4.3. SZ angular power spectra 

This section focuses on the frequency-dependent SZ an- 
gular power spectrum, mainly near 150 GHz where there 
are recent measurements by ACT (.F owler e t al.l |2010[) 
and SPT (|Lueker et al.ll2009l : IHall et al.ll2009lV We first 
compare the models below, then discuss how they scale 
with (78 in Section [4.4l and finally discuss the constraints 
from recent observations in Section [S] 

Figure [S] shows our four models and the KS02 model 
at 148 GHz, which is one of three central observing fre- 
quencies of ACT. Relative to the standard model at 
/ = 3000, the adiabatic model is 41% larger, the low- 
/gas model is 27% smaller, and the nonthermal20 model 
is 37% smaller. The trend in the model amplitudes is 
the same as found for the TSZ effect, since it contributes 
several times more power than the KSZ effect at this fre- 
quency. In comparison, the KS02 template now has less 
power than the standard model at / > 2000 because it is 
based only on the TSZ effect. However, the differences 
between the two models are < 10% for the multipole 
range of current interest (/ w 3000 ± 1000), where the SZ 
power begins to exceed that of the CMB but not yet over- 
whelmed by the rising power from point sources. Thus, 



the standard and KS02 templates should give similar re- 
sults when used to fit for the SZ contribution to the CMB 
angular power spectrum near this frequency. 

4.4. Dependence on erg 

In order to place constraints on erg, we first 
need to know how the TSZ, KSZ, and SZ angular 
power spectra scale with this cosmological parameter. 
iKomatsu & Seljak (2002) suggested that the TSZ power 
scales approximately as Ctsz oc Ug. Using perturbation 
theory, iVishniad (|1987D showed that the large-scale OV 
angular power spectrum goes as Cqv oc (t|. With the ad- 
dition of nonlinear contributions, the KSZ power should 
have a stronger dependence. To quantify how the SZ 
templates scale with erg, we define the (q = SZ, TSZ, 
KSZ) scaling amplitude. 



Cqa,ag =0.8) 



V0.8/ 



(9) 



which is assumed to be a powerlaw relation with dg. In 
general, the scaling amplitude and the scaling index. 



din 
dins 



(10) 



where s = (7g/0.8, are also dependent on multipole Z, 
frequency v, and astrophysical and cosmological param- 
eters. 

We can rescale the four models to arbitrary as with- 
out having to run additional simulations. Each model 
is constructed from gas in identifiable halos and in the 
IGM, and their contributions are rescaled as follows. 
At / ^ 3000, gas in identifiable halos (with MpoF > 
7 X 10^2 H-^Mq at z < 3) contribute ~ 95% and ~ 30% 
to the TSZ and KSZ power spectra, respectively. Since 
the power from this component is domi nated by the one- 
halo term (jKomatsu Sz Kitavamall9"99| ) and proportional 
to the halo mass function, we can rescale it with the mul- 
tiplicative factor. 



»^q(^cr8) 



Y,CqAM,z,l,as = 0.8)u;,(M,z,crg)u^,(crg) 

i 



which is calculated by weighting the contribution of each 
halo i by the ratio of the number density per unit mass. 



w{M,z,as) 



dn 



dn 
dM 



(12) 



(M, Z, erg ^ 0. 



We use the IJenkins et alj ()20Q1D fitting formula for the 
halo mass function, which agrees well with the simulated 
halo abundance. The velocity scaling factor usz for the 
SZ power is a linear combination of utsz = 1 and mksz = 
crg/0.8. Gas in the IGM at z < 3 and all gas at 3 < 
z < 10 together contribute ~ 5% to the TSZ power, but 
~ 70% of the KSZ power at I ~ 3000. For this low- 
density component, we assume the linear- regime scaling 
factors ttsz = (cg/O.S)^ and rxsz = (o-g/0.8)''. Since 
the power is dominantly from the OV effect, rsz ~ t'ksz 
at all frequencies. 
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(with Q!L = 4) and a nonlinear term Cnl (with onl > 4), 
the kinetic scaling index can be written as 



Fig. 9. — Top: Scaling index axsz for the TSZ angular power 
spectrum for the standard (black solid), adiabatic (red long-dash), 
low- /gas (blue short-dash), nonthermal20 (green dotted), and KS02 
(orange short-long-dash) models at I = 3000. Middle: Scal- 
ing index oksz for the KSZ angular power spectrum. Bottom: 
Frequency-dependent scaling index asz at 148 GHz. 

Figure ini shows the effective scaling index as a function 
of as 'At I = 3000. For the TSZ component, we find 
7 ^ Q!TSZ ^ 9 for 1 > CTg > 0.6; at erg ~ 0.8, axsz ~ 7.7 
(adiabatic), 8.1 (standard), 8.2 (nonthermal20), and 8.5 
(low- /gas)- The thermal scaling index increases in models 
where massive or higher redshift halos make a relatively 
larger contribution to the power. Considering Ctsz as 
the sum of halo terms Chaio.i and an IGM term Cigm, 
the thermal scaling index can be written as 



din (a;iGMS^ + ^ 



axsz 



d In s 



(13) 



where Xhaio.i = C'haio,j(s = 1)/Ctsz(s = 1) and xigm = 
Cigm(s = 1)/Ctsz(s = 1). The halo abundance scaling 
w{M, z,as), defined in Equation (jl2p . rapidly increases 
(decreases) with mass and redshift at fixed erg greater 
(less) than 0.8. As discussed in Section I4TT1 both mas- 
sive halos (M500 ^ 3 X 10^^ Mq) and high redshift ha- 
los contribute increasingly larger fractions of the thermal 
power going from the adiabatic model to the low- /gas 
model. This explains why axsz is smallest for the adia- 
batic model and largest for the low- /gas model. At any 
given mass and redshift, w monotonically increases with 
as, but its slope, dw/das, gradually decreases, resulting 
in the generally observed relation daxsz/dcrg < 0. In all 
of the models, axsz at as ~ 0.8 decreases by approxi- 
mately 0.5 going from I — 3000 to I = 8000 since low 
mass halos become relatively more important at smaller 
angular scales. 

For the KSZ component, we find axsz ~ 4.8 (low- /gas), 
4.8 (nonthermal20), 4.9 (standard), and 5.0 (adiabatic) 
at as ~ 0.8. This trend is different from that for the TSZ 
effect. Considering Cksz as the sum of a linear term Cl 



ttKSZ 
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where Xh = Cl(s = 1)/Cksz(s = 1) and xnl = Cnl(s = 
1)/Cksz(s = 1)- Since all of the models share the same 
linear component, aKSZ only depends on the nonlinear 
variables. On one hand, the nonlinear power ratio xnl 
increases in models with larger gas fractions, being small- 
est for the low- /gas model and largest for the adiabatic 
model, as seen in Figure [71 On the other hand, the non- 
linear scaling index onl has the opposing trend, being 
lowest for the adiabatic model and highest for the low- 
/gas model, as with axsz in Figure [HI In our models, 
the trend in axsz is explained by the fact that the dif- 
ferences in xnl slightly outweigh the differences in aNL- 
At I — 3000, there is weak variation with erg, since the 
linear component is generally more dominant and ol is 
independent of ug (assuming homogeneous reionization) . 
At / ^ 8000, the larger nonlinear ratio results in axsz 
increasing by approximately 0.3 at erg ~ 0.8, and there 
is larger variation with erg. 

Figure [5] also shows asz at 148 GHz. Since axsz and 
axsz are different functions of erg , asz will be frequency- 
dependent. In the nonrelativistic limit, the SZ power is 
given by Csz('^) — C'TSz(t^) + C'ksz, and the SZ scaling 
factor and index can be written as 



^sz = a;xsz^TSz + a;Ksz^Ksz, 



(15) 



din ixTSZS°"^^^ + XKSZS°"^^^) 

asz = , (16) 

d m s 

respectively, where xtsz{i^) = C'xsz(s = 1)/C'sz(s ~ 1) 
and XKszH = Cksz(s = 1)/Csz(s = 1)- At erg w 0.8, 
we find asz ~ 7.3 (nonthermal20), 7.4 (adiabatic), 7.5 
(standard), and 7.6 (low- /gas)- The trends for asz are 
generally similar to those for axsz, because the thermal 
component contributes several times more power than 
the kinetic at this frequency. However, the nonthermal20 
model is an exception. Even though it has the second 
largest axsz amongst the models, it has the lowest asz 
because it also has the lowest thermal power ratio xxsz- 
The SZ scaling index does not vary monotonically with 
erg, and the relation can be complicated in general- To- 
ward lower erg, axsz gradually increases, but the lower 
aKSZ increases in weight. Toward higher erg, axsz grad- 
ually decreases, but it also has more weight. 

We also calculate axsz for the KS02 template by 
repeating the ha l o mo del calculation described in 
iKomatsu fc SeliakI ([200I - At erg « 0-8 and I = 3000, 
axsz ~ 7.9 is higher than the previously suggested value 
of approximately 7. With the addition of the KSZ com- 
ponent, asz should be similar to that from the standard 
and adiabatic models. 

5. DISCUSSION 

ACT and SPT have used our standard model to fit for 
the SZ contribution to the CMB angular power spectrum 
on arcminute scales. Figure |5] shows the constraints on erg 
based on the measured ^sz,std, which we convert to an 
SZ power at / = 3000 and 148 GHz to compare with the 
models. The following comparison is based on published 
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Fig. 10. — Constraints on erg from the SZ power at Z = 3000 
and 148 GHz, assuming the standard (black solid), adiabatic (red 
long-dash), low- /gas (blue short-dash), nonthermal20 (green dot- 
ted), and KS02 (orange short-long-dash) models. The ACT (95% 
CL) upper l imit (Fowler ot al. 2010) and the SPT best-fit value 
iLueker et al.,.2009 . Icr uncertainty in gray) for the SZ power are 
derived using the reported Agz gxD- 

SZ data points and not done rigorously by fitting the 
obser ved CMB power spectrum. 

Fowler et alj (|2010( ) place an upper limit (95% confi- 
dence level) of Asz,STD < 1-63 based on ACT obser- 
vations at 148 GHz. This corresponds to < 16.8 /iK^ 
at / = 3000 and 148 GHz. We find as < 0.85 (stan- 
dard), 0.82 (adiabatic), 0.89 (low- /gas), and 0.91 (non- 
thermal20). Since the adiabatic model is known to be in- 
correct, it is safe to omit the constraints from this model. 
These bounds remain consistent with t he most recent 
value from the WMAP 7-year results (jKomatsu et al.l 
|201C1[) . but they are not highly constraining. With only 
one of three observing frequencies used in this first anal- 
ysis, the degeneracy between the SZ and point source 
contributions to the CMB angular power spectrum can 
not be broken to obtain best-fit values. Upcoming ACT 
measurements, which will also include observations at 
approximately 220 and 280 GHz, will provide much 
stronger c onstr ai nts. 

Luc ker et al.l ()2009t ) report a best-fit value of 
Asz,STD = 0.42 ± 0.21 (at 153 GHz) for SPT obser- 
vations near 150 and 220 GHz. This corresponds to 
4.3 ± 2.2 ^K^ at / = 3000 when converted to 148 
GHz. We find erg = O.TligjJe (standard), 0.68^^:^6 
(adiabatic), 0.74j:°:g^ (low-/gas), and 0.76tJ^!]^ (non- 
thermal20). The central values from this first analy- 
sis are general ly low compared to the WMAP 7-year 
measurements (jKomatsu et al .11201 Of ), but still within la 
for the low- /gas and nonther mal20 models and w ithin 
2(7 for the standard model. iLueker et al.l (|2009f) find 
o-g = 0.773 ± 0.025 in a joint analysis w ith the WMAP 
5-year constraints ([Dunklev et al.ll2009[ ). We expect a 
higher as, closer to 0.8, if the same analysis were repeated 



assuming the nonthermal20 or low- /gas models instead of 
the standard model. Upcoming SPT analysis with more 
survey area and with a third observing frequency of 90 
GHz, will provide more conclusive results. 

In addition, ACT d Hincks et al.' "200?) and SPT 
()Staniszewski et al.l l2009i .Plagge et al. 2010.) have de- 
tected the SZ efi'ec t directly for mass ive galaxy clus- 
ters. In particular, iPlagge et al.l ()201Cl[ ) found that the 
stacked radial profile from 15 clusters has a very sim- 
ilar shape to tha t predict e d by the best-fit pressure 
profile from Arna ud et al.l ()2009| ). Using the best- fit 
GNFW parameters provided for the individual clusters, 
we compared their angular profiles against a fair sub- 
set (M500 > 5 X 10" Mq and z < 0.4) of cluster profiles 
from the standard model, finding good agreement in both 
shape and amplitude. If we correct the X-ray masses for 
the assumed nonthermal hydrostatic support, then the 
SPT cluster profiles are also in good agreement with the 
nonthermal 20 model. 

Recently, iKomatsu et"all ()2010[ ) analyzed the WMAP 
7-year maps and report a significant (8a) detection of 
the SZ effect at the locations of known galaxy clus- 
ters from the REFLEX catalog (Bdhringcr ct al. 2004). 
They stacked 742 clusters with z < 0.45 and calculated 
the average projected SZ angular profile from the CMB 
maps. The measured signal is found to be approximately 
40% smaller than expec tation based o n the best-fit pres- 
sure profile from Arnau d~et al.l (|2009( ). This discrepancy 
could be due to the possibility that the GNFW pres- 
sure profile derived jointly from X-ray observations and 
numerical simulations could be overestimating the pres- 
sure at R/R500 ^ 1 or R/R500 ^ O-l- Unfortunately, 
we can not fairly compare our models since the stacked, 
weighted signal depends on detailed properties of the 
WMAP beam and noise, which are unknown to us. How- 
ever, the discrepancy may not originate with the adopted 
pressure profile, but rather with the highly uncertain 
scahng relations (e.g. i?5oo — Tx) used to calculate the 
expected signal. 

6. CONCLUSION 

We have modeled the SZ effect by postprocessing a 
simulation of the large-scale-structure of the universe to 
include gas in dark matter halos and in the filamentary 
IGM. All free electrons are accounted for up to the epoch 
of cosmic reionization, which is assumed to occur homo- 
geneously at z = 10 (tt ~ 0.08). See Schg al et al.l ([20T0I ) 
for additional details and discussion. 

We calculate the first-order TSZ (Compton y) and 
KSZ (Doppler b) signals, and the frequency-dependent 
SZ temperature fluctuations. In the latter, the ther- 
mal and kinetic components have higher-order relativis- 
tic corrections and are nonlinearly coupled. Sky maps 
with pixel resolution of 0.4 arcminutes are made by trac- 
ing through the simulated electron distribution. 

Four models for the halo gas distribution are compared 
to study how gas physics influence the SZ temperature 
fluctuations and angular power spect ra: adiabati c, stan - 
dard, low- /gas, and nonthermal20. iBode et al.l (|2009f ) 
previously calibrated the standard model to reproduce 
the stellar and gas fractions and X-ray scaling relations 
measured from low redshift clusters and groups. Relative 
to the standard model, the other models differ by approx- 
imately ±50% (TSZ), ±20% (KSZ for Zreion = 10), and 
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±40% (SZ at 148 GHz) at / = 3000 and erg = 0.8. 

Wc also calculate the dependence of the angular power 
spectra on as- The templates for (Tg — 0.8 are rescaled 
to arbitrary ag without having to run additional simula- 
tions. To quantify the dependence of C; on ag, a pow- 
erlaw scaling relation Ci oc (as/OS)" is assumed, where 
a = a{l, (Tg) is the effective powerlaw scaling index. For 
the collection of considered models, we find 7 < axsz ^ 
9, 4.5 < axsz < 5.5, and 6.5 < asz(148 GHz) < 8 for 
0.6 < (78 < 1 at / = 3000. Care should be taken when 
applying the reported scaling indices to other multipoles, 
frequencies, or to other gas models. 

We summarize some interesting features of the SZ an- 
gular power spectrum in our models: 

1. The SZ angular power spectrum decreases in am- 
plitude as gas mass and binding energy is removed 
through star formation, and as gas is pushed out 
to larger radii by energetic feedback. 

2. Nonthermal pressure support can have comparable 
effects as energetic feedback on lowering the ther- 
mal contribution to the SZ power (for reasonable 
choices of parameters). 

3. The broad peak in the power spectrum shifts to 
larger angular scales as massive or low redshift ha- 
los make a relatively larger contribution to the total 
power. 

4. The TSZ scaling index axszCcs) increases as 
massive or high redshift halos make a relatively 



larger contribution to the power. It also decreases 
with both I and erg. 

5. The KSZ scaling index Q!ksz(^,CT8) increases with 
the nonlinear power ratio Cnl/Cksz- It increases 
with I, but has only weak variation with as for 
homogeneous reionization. 

6. The frequency-dependent SZ scaling index 
oisz{l,crs,i^) increases with the thermal power 
ratio Ctsz/C'sz- It does not have to vary mono- 
tonically with as and can be complicated in 

general. 
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